Exponentially Fragile PT-Symmetry in Lattices with Localized Eigenmodes 
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We study the effect of localized modes in lattices of size N with parity-time (PT) symmetry. 
Such modes are arranged in pairs of quasi-degenerate levels with splitting S ~ exp _JV/ ^ where £ is 
their localization length. The level "evolution" with respect to the PT breaking parameter 7 shows 
a cascade of bifurcations during which a pair of real levels becomes complex. The spontaneous PT 
symmetry breaking occurs at "/j>t ~ min{<5}, thus resulting in an exponentially narrow exact VT 
phase. As N/£ decreases, it becomes more robust with 7-pr ~ 1/N 2 and the distribution P(7pt) 
changes from log-normal to semi-Gaussian. Our theory can be tested in the frame of optical lattices. 



Introduction - Parity (P) and time -reversal (T) 
symmetries, as well as their breaking, belong to the 
most basic notions in physics. Recently there has been 
much interest in systems which do not obey P and T- 
symmetries separately but do exhibit a combined PT- 
symmetry. Examples of such PT-symmetric systems 
range from quantum field theories to solid state and 
classical optics [IJ|3|31II1|5J|6J|71IE1I5]- A "PT- 
symmetric system can be realized in optics, by creating 
a medium with alternating regions of gain and loss, such 
that the (complex) refraction index satisfies the condi- 
tion n*( — x) — n(x) [SI [HI IS] ■ This condition implies 
that creation and absorption of photons occur in a bal- 
anced manner, so that the net loss or gain is zero. Such 
synthetic PT metamaterials show unique characteristics 
such as "double refraction" and non-reciprocal diffraction 
patterns, which may allow their use as a new generation 
of unidirectional optical couplers or left-right sensors of 
propagating light [BJ. In the paraxial approximation the 
classical wave equations reduces to a Schrodinger equa- 
tion with a fictitious time, related to the propagation 
distance, and with the refraction index playing the role 
of the potential. We use below the terminology of the 
Schrodinger equation, while keeping in mind applications 
to optical systems. 

A PT-symmetric system can be described by a phe- 
nomenological " Hamiltonian" Ti. Such Hamiltonians 
may have a unitary time evolution and a real energy 
spectrum, although in general arc non-hermitian. Fur- 
thermore, as some parameter of Ti changes, a sponta- 
neous VT symmetry breaking occurs, at which point the 
eigenfunctions of Ti cease to be eigenfunctions of the PT- 
operator, despite the fact that Ti. and the PT-operator 
commute pQ. This happens because the PT-operator is 
not linear, and thus the eigenstates of Ti. may or may 
not be eigenstates of VT. As a consequence, in the bro- 
ken VT -symmetry phase the spectrum becomes partially 
or completely complex. The other limiting case where 
both Ti. and VT share the same set of eigenvectors, cor- 
responds to the so-called exact VT -symmetric phase and 
the spectrum is real. 



In this Letter we investigate the spontaneous PT- 
symmetry breaking scenario in a wide class of systems 
supporting localized states. Such states are ubiquitous 
in macroscopic systems. They can reside on impuri- 
ties or at the edges of an otherwise perfect lattice of 
finite size. Therefore, in order to understand the PT- 
symmetric phase of a macroscopic system, it is impera- 
tive to consider localized states. At the same time, we 
note that even fifty years after the seminal work of An- 
derson 10 . localization continues to be a thriving area of 
research, not only for solid-state physics, but also to other 
fields including ultra-cold atoms, acoustics, microwaves 
and classical optics. We therefore expect that our study 
linking the newly developed area of VT materials with 
the field of localization will contribute to understanding 
fundamental aspects of modern physics. 

Localization is particularly pronounced in one- 
dimensional (ID) systems and has been studied exten- 
sively in the past [TTJ. We show that in case of PT- 
symmetric lattices of size N which can support localized 
modes due to disorder or impurities or even due to bound- 
aries (surface states), the mechanism that triggers the 
transition from real to complex spectrum is level cross- 
ing between a pair of modes having the smallest energy 
spacing. Due to the P-symmetry, this pair of states has 
a double hump shape and the energy splitting between 
them is Si ~ exp(— Zo/£) where £ is the localization length 
and Iq is the distance between the two humps (for disor- 
dered lattices Iq ~ N) . We find that the value of the VT 
symmetry breaking parameter 7 at the transition point is 
j-pt ~ Si , thus indicating that the exact VT -symmetric 
phase is exponentially small in the limit Iq/^ 3> 1. In 
contrast, for Zq/£ "C 1, we find that the smaller level 
spacing scales as A m j n ^ 1/N 2 . This is also reflected in 
the distribution V(j-pt) which changes from a log-normal 
towards a semi-Gaussian as N/^ decreases. 

Two PT-symmetric Impurities - It is instruc- 
tive to start with the simple example of a pair of PT- 
symmetric impurities implanted into an otherwise perfect 
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infinite lattice. The system is described by the equation 

- l/'n+l - lpn-1 = (E - e n )lp n , (1) 

where ip n is the eigenfunction amplitude at site n, e n = 
for n ^ ±1, and e±i = —(3 ± 17, with (3 and 7 being real 
and positive. We are looking for the bound states: 

( Ae kn , n < -I 

i> n = < Be kn + Ce~ kn , -l<n<l (2) 
[ De- kn , n > I 

with Re [A;] > and E = —2 cosh A;. Matching the wave 
function at the sites n = ±Z, one obtains four equations 
for the amplitudes A, B, C, D. Equating the determinant 
to zero yields the transcendental equation for k: 

sinhfc=^±i[- 7 2 + (/3 2 + 7 2 )e- 4fc f /2 . (3) 

For 7 = and (31 3> 1, one finds two bound states with 
energies E± — E ±^6i, where E = — \J\ + (3 2 is the en- 
ergy of a localized state on a single isolated impurity, and 
61 = 6(1) = (2(3 2 / |Eb|)e - ^ is the exponentially small 
energy splitting term for the two-impurity problem. The 
point we want to emphasize is that for (31 3> 1 already 
an exponentially small 7 leads to complex values of k 
and E, thus, breaking the PT-symmetry. The mecha- 
nism for this breaking is level crossing: as follows from 
Eq. ([3]), when 7 reaches the value 7 -pr ~ /3e~ ,3/ , the 
two (real) eigenvalues become degenerate. For 7 > 7 -pr 
they branch out into the complex plane, displaying near 
the branch point the characteristic behavior Im[£'±] oc 
ib^/7 2 — 7-pT-. This square root singularity seems to be 
a generic feature of the PT -symmetry breaking. 

The eigenfunctions of the above Hamiltonian also un- 
dergo characteristic changes as 7 increases. A finite 7 
breaks the P-symmetry of the Hamiltonian but, as long 
as 7 < 77377 the PT-symmetry of the eigenfunction is 
preserved, so that ^* = ±-0_„. This implies that the 
"dipole moment", D — J2 n =-N n IVvJ 2 , of an eigenstate 
is zero. For 7 > 7 -pr the eigenstates acquire a finite 
dipole moment. Below we shall use D as one of the sig- 
natures of the PT-symmetry breaking. 

The disordered PT-symmetric Chain - We con- 
sider next a ID disordered PT-symmetric chain and 
demonstrate that under a disorder increase, the PT- 
symmetric phase is gradually destroyed. For sufficiently 
strong disorder this phase shrinks to an exponentially 
narrow region, even for a comparatively small system. 

The chain is described by the Hamiltonian of Eq. 
0, where now all e n are random complex numbers, 
£n = /3n + *7«, with the constraint e„ = e*_ n . One can en- 
visage various possibilities for randomness, either in (3 n or 
7„, or both. Below we present results for the case where 
(3 n (for n > 0) are uniformly distributed on the interval 
[—(3/2; (3/2] and 7 „ = 7 = const, (for n > 1). It is cru- 
cial, for the PT-symmetry, to implement the constraint 
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FIG. 1: Pairs of exponentially localized modes in a ID-chain 
with PT-symmetric disorder (a-c) and surface states in a 
PT-symmetric periodic chain (d) yielding the smallest energy 
splittings Si for (3 — 3. For 7 < 7-pr (b) the two eigenfunc- 
tions (blue and red; the imaginary part is shown in the in- 
set) are complex but PT symmetric and their absolute values 
remains equal and symmetric [coinciding on the black line 
in (a)]. For 7 > j-pr (c) the eigenfunctions (magenta and 
cyan) are no longer PT-symmetric and shift their weight to- 
wards separate sides of the chain (a), (d) Surface states (blue 
and green) showing exponential localization. The red dashed 
[solid] lines in (a) [(d)] are guides to the eye, pointing out the 
exponential localization. 



(3 n = (3- n and 7„ = — 7- n (the latter implies 70 = 0). 
This constraint introduces a peculiar long-range correla- 
tion. To clarify the picture we start with the Hcrmitian 
limit 7 = 0, and assume a long chain, such that even 
eigenstates in the band center are localized, i.e. their lo- 
calization length is smaller than the system size (27V + 1). 
Imagine for a moment that the chain is cut in half, at its 
center n = 0. Then a typical state in one half of the 
chain would be localized, with some localization length 
£, far away from n — 0, say, near site i > 1. This state 
has its counterpart in the other half of the chain, near 
site —I. In the full, connected chain this pair of states 
has an exponentially small overlap, at the center of the 
chain, yielding two real eigenstates of the entire chain. 
Each of these eigenstates (one symmetric, the other an- 
tisymmetric) has two peaks, near the sites n = ±Z. The 
energy splitting between the two eigenvalues is exponen- 
tially small, 6(1) ~ e^ 21 ^, in complete analogy with the 
example of the two impurities in a perfect chain. 

Thus, the eigenstates in a P-symmetric disordered 
chain are organized into pairs (doublets) of energy split- 
ting 61 < 6 2 < ■ ■ ■ ■ The energy splitting between the 
two states of a doublet is exponentially small, while the 
energy separation between consecutive doublets is much 
larger, of the order of level spacing, A, in half chain. The 
pair associated with the minimum splitting, 61 is likely to 
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originate from states localized far away from the origin 
of the chain (n = 0) and with energies at the band edges 
(small £) . As 7 is switched on the eigenstates of each pair 
will initially preserve their PT-symmetric structure (see 
Fig. la,b). At 7 = -f-pr — Si, the two levels associated 
with Si will cross, breaking the PT-symmetry (see Fig. 
2a). As 7 > 7-pr these modes cease to be eigenstates of 
the PT-operator. Instead, the weight of each of them is 
gradually shifted towards one of the localization centers. 
An example of such pair associated with 5\ is reported in 
Fig. la,c. For larger 7 the next doublet, with splitting 
82 will come into play, creating a second pair of complex 
eigenvalues for 7 ~ 82 (see Fig. 2a), etc. 

The above qualitative considerations apply to a single, 
realization of the random potential. A full theory must 
be formulated in statistical terms and deal with proba- 
bility distributions. For instance, the critical value j-pt 
at which the VT -symmetry is broken, fluctuates from 
one realization to another and the appropriate question 
pertains to the distribution P(j-pr)- As has been argued 
above (see inset of Fig. 2a), in the strong disorder regime 
7737- 8%, and thus the problem reduces to the study of 
the distribution V{5\). There are several sources of fluc- 
tuations in 81: fluctuations in the position and energy 
of the relevant localized states, as well as what can be 
termed "fluctuations in the wave functions" . By this we 
mean that a localized wave function exhibits strong, log- 
normal fluctuations in its "tails" , i.e. sufficiently far from 
its localization center [12] . This latter source of fluctua- 
tions appears to be the dominant one and it immediately 
yields a log-normal distribution for Si (see Fig. 2b), since 
81 is proportional to the overlap integral between a pair 
of widely separated and strongly localized states. 

The aforementioned scenario of bifurcations, i.e. of 
the consecutive branching of pairs of eigenvalues into the 
complex plane, can also be made more quantitative. Con- 
sider, the separation (on the 7-axis) between the first bi- 
furcation (7 = 7-pr) and the next one. This separation, 
8-f, is proportional to (8 2 — 81). Assuming that local- 
ized states, associated with the doublets, are randomly 
located, and ignoring for the moment fluctuations in the 
energy of these states, one immediately obtains a Poisson 
distribution for s = log 82 — log 8% , with the average spac- 
ing 2/£ on the log 7 scale between the bifurcation points 
(see inset of Fig. 2b). This result, with £ being replaced 
by an appropriate average, remains valid also when we 
account for the energy fluctuations. 

Our considerations can be extended to the N/£ <C 1 
limit, when the states responsible for 81 are extended 
over the entire chain. In this case the picture of doublets 
with exponentially small energy splittings is not valid 
and 7-pr becomes of the order of the minimal level spac- 
ing, A m j n , in the corresponding Hcrmitian problem. This 
statement follows from perturbation theory, with respect 
to 7. The unperturbed (i.e. 7 = 0) energy levels are 
real, and are separated by intervals of order 1/N 2 at the 
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FIG. 2: (a) Bifurcations for the dipole moment D and for 
the imaginary part Im_E of energy levels for a PT-symmetric 
ID chain with 7 =const. and /3 n given by a box distribution 
[— §; § ] for N/£ 3> 1. The first two bifurcations (correspond- 
ing to splittings 81 and 82) are shown. The square-root behav- 
ior at the branching point (see text) is indicated with dashed 
cyan lines. Inset: a linear relation between 81 and 7pt is 
evident for almost 10 orders of magnitude, (b) Distribution 
V(x) of x = (log^i — (log5i))/<7 (er is the standard devia- 
tion) for various disorder strengths /3. In the limit of large j3 
the distribution becomes log-normal. Inset: The distribution 
V(s) of s = log(<52) — log(<5i) reported in a semi-logarithmic 
plot. A Poisson distribution is approached as f3 increases. 



band-edges (at the center of the band the separation is of 
order 1/N), so that A m ; n ~ 1/N 2 . Finite 7 leads to level 
shifts proportional to 7 2 (the first order correction van- 
ishes due to PT-symmetry) and for 7 = j-p-r ~ A m i n the 
perturbation theory breaks down, signaling level crossing 
and the appearance of the first pair of complex eigenval- 
ues. Thus, the energy scale for the PT-threshold in the 
N/£ < 1 limit (j VT ~ 1/iV 2 ) widely differs from that 
for N/l; 3> 1 (civr — e~ 2N /^). However, the "bifurca- 
tion scenario" , with characteristic square-root branches 
for the complex eigenvalues, holds in both cases (again, 
with the completely different energy scale for the inter- 
vals between the consecutive bifurcations). Our numeri- 
cal results presented in Fig. 3a are in perfect agreement 
with these considerations. 

We study now the distribution P{^-pt) in the limit 
iV/£ <C 1. We invoke perturbation theory with respect to 
the perfect lattice. The perturbative scenario, indicates 
that weak disorder will cause a small shift of the energy 
levels. Thus the new level spacing becomes A m ; n ± AS, 
where the sign + (— ) refers to the upper (lower) band- 
edge. Regardless of the sign of Ai5, the minimal level 
spacing (in first order perturbation theory) is 
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with the coefficients A„ = sin 
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If f3 n were Gaussian distributed, it would be immediately 
clear that the distribution of 8, P{8), is a semi-Gaussian 
for 5 < A m ; n . This should remain approximately true 
also for the box-distribution, employed in our numerics, 
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FIG. 3: (a) Same as in Fig. 2a but now for "weak" disor- 
der N/£ <C 1. The same bifurcation scenario is observed, 
(b) The distribution P(A) of the rescaled minimal energy 
split A = AS /a (a is the standard deviation), for various 
disordered strengths, all of them being in the weak localized 
regime. The distribution has an upper cutoff at A = 0. A 
Gaussian distribution is shown also by a red dashed line. In- 
set: The same data in a semi-logarithmic manner vs. A 2 . 



if the number of terms in the sum is sufficiently large. 
Fig 3b confirms this expectation. 

Periodic PT-symmetric Potentials - Let us 
briefly address the problem of PT-symmetry breaking 
for a periodic complex potential. This question has been 
raised in |6j, for a potential V(x) — 4(cos 2 £-MVb sin 2a;), 
where it was argued that the critical value, of the PT- 
threshold was Vq = 1/2. The presence of the real part, 
A cos 2 a; (A =^ 0), is crucial for this result. Indeed, it 
was pointed out in [6J that a purely imaginary periodic 
potential, V(x) = zVb sin 2Ar+1 (x), treated in Ref [2], has 
"zero PT-threshold" , i.e. it cannot have an entirely real 
spectrum. Another example of a periodic potential with 
zero PT-threshold was provided in Ref [3]. 

In a periodic system of finite extent one usually en- 
counters localized states (surface states) at the bound- 
aries of the sample jl3J . We have found that these states, 
which were not addressed in Ref. [B], are crucial for the 
correct evaluation of the PT-threshold in a finite sys- 
tem. We illustrate the importance of the surface states 
by taking the example of a tight binding periodic poten- 
tial, with three sites per unit cell. The Hamiltonian is 
still that of Eq.Q, but with on-site energies: £ n= 3/ = 
0; e„=3;±i = ± ry where I = 0, ±1, ±2, . . . , ±N/3. For 
7 = (and N — > oo) the spectrum displays two en- 
ergy gaps, whose width (for (3 <§; 1) is 2/3/3. The ex- 
istence of gaps suggests, in analogy with 6 , that the 
PT-symmetric phase in this model will exhibit some ro- 
bustness, as long as 7 is small (7 <§; /3). This expectation, 
however, is not born out due to the surface states. For 
7 = there is a pair of surface states, exponentially de- 
caying away from the sites ±7V. For large but finite N 
these surface states have a small overlap near the cen- 
ter of the chain n = 0, yielding a doublet with an ex- 
ponentially small energy splitting. In complete analogy 
with the two-impurity problem, already an exponentially 



small 7 ~ e~^ N will therefore break the PT-symmetry. 
This example shows that the N — > 00 limit can be quite 
subtle, as far as the PT-threshold is concerned. If one 
starts directly with N = 00, one obtains a finite PT- 
threshold, 7-pr — /3- If one starts, however, with a large 
but finite N and then takes the A~ — > 00 limit (which 
is the correct physical limit), then one ends up with 
Ivt = 0, due to the existence of surface states. 

Conclusion - We have studied a ID PT-symmetric 
chain with disorder. The PT- symmetric phase turns out 
to be very fragile. For a sufficiently long chain, this phase 
exists only for exponentially small values of the imaginary 
part of the potential j-pj- ~ e~~ N ^ beyond which the PT- 
symmetry is broken (here N and £ are the system size 
and the localization length respectively) . Our model can 
be extended in various directions. For instance, we have 
checked that our main results do not change if random- 
ness is introduced in both the real and the imaginary 
parts of the potential. We have also briefly discussed the 
periodic PT-potential and pointed out the importance of 
surface states in breaking PT-symmetry. Our main con- 
clusion is within a generic ID system which supports lo- 
calized modes, the threshold for PT-symmetry breaking 
is exponentially approaching zero with increasing size. 
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